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We generalize the calculation of Casimir interaction between topological insulators with opposite 
topological magnetoelectric polarizabilities and finite surface band gaps to finite Temperature cases. 
We find that finite temperature quantitatively depress the repulsive peak and enlarge the critical 
surface gap rric for repulsive Casimir force. However the universal property rrica ^ 1/2 is still valid 
for various oscillation strength, temperature region and topological magnetoelectric polarizabilities. 
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I. INTRODUCTION 

Exploration for exotic physical properties about topo- 
logical protected quantum states is an important theme 
of current condensed matter physics. The recently dis- 
covered topological insulator (TI) is such a quantum 
state. The three dimensional topological insulator has 
a bulk gap like an normal insulator, however the sur- 
face state of this material is gapless^^^^ , and such a 
gapless spectrum together with odd Dirac cones on TI 
surface are topological protected by the time-reversal 
symmetry^ There are many interesting phenome- 

nas (predictions) related to this novel material, such as 
the topological magnetoelectric effect , electric charge 
induced magnetic monopole^^'^^, optical Kerr and Fara- 
day rotation , surface 1/2 quantum Hall effect^^'^^, 
et.al. 

Casimir effect is a quantum effect arising from zero- 
point energy fluctuation of vacuum, the seminal work of 
H. B. G. Casimir found that two parallel uncharged 
metallic planes will emergence an attractive force. Be- 
fore investigated in TI systems, casimir force has been 
proposed to be repulsive for some special conditions. 
For instance, it is proposed that Casimir force is repul- 
sive if special geometry has been considered ^ , it is also 
reported that Casimir interaction between metamateri- 
als maybe repulsive'^^'^'^, experimental evidence shown 
that high-refractive liquid^'^ between dielectrics will in- 
duce repulsive Casimir force. 

Recently, A. G. Grushin and A. Cortijo proposed^^'^'' 
that Casimir interaction between TIs with opposite topo- 
logical magnetoelectric polarizability is repulsive while 
the distance between TIs tends to zero. Their analyza- 
tion is based on the topological quantum field descrip- 
tion of TI , Stopo = Q^/(47r)^ / (fxdtOE • B, where a is 
the fine structure constant, 6 = (2n + l)7r is the topo- 
logical magnetoelectric polarizability, E and B are elec- 
tric and magnetic field respectively. Such a topological 
quantum field description is exact, however, the repul- 
sive Casimir force will be suppressed by conducting sur- 
face fermions. In order to deduce the influence of surface 
fermions, one need to open a surface band gap by adding 



a magnetic coating on TI. We analyzed the Casimir in- 
teraction between TIs for finite surface band gap at zero 
temperature . We found that a critical surface band 
gap rUc is essential for repulsive Casimir interaction, and 
such a surface band gap can be estimated by rrica ~ 1/2, 
where a is the distance between TIs. 

For practical measurement, the effect of temperature 
is always need to be considered. In this paper, we calcu- 
lated the Casimir interaction between TIs with opposite 
topological magnetoelectric polarizability at finite tem- 
perature, we found that the general relation rrica ~ 1/2 
is still valid. 

This paper is organized as follows: In Sec. II, we derive 
an effective action of surface electromagnetic field by 
integration out the contribution of surface fermions with 
finite surface band gap. From the effective action, we 
deduce the Maxwell equations with boundary corrections 
and Fersnel coefficient matrix in Sec. III. In Sec. IV, we 
calculate the Casimir energy between TIs by Lifshtz 
formula, then we present the scope of repulsive Casimir 
force for different temperatures, surface band gap and 
topological magnetoelectric polarizabilities. Conclusions 
are given in Sec.V. 



II. EFFECTIVE ACTION AT FINITE 
TEMPERATURE 

Let us formulate the model, in the vacuum and bulk 
of TIs, the action of electromagnetic field can be written 
as: 

s = -llm'.{.E'-lB'). (1) 

where E and B are electric and magnetic field, e and fi 
are permittivity and permeability of TI in the bulk and 
equal to 1 in the vacuum. 

The topological nontrivial term a/ (47r)^ / d^xdtOE • B 
can be modeled by massive surface Dirac fermions: 

Sd= I d^x^ [ij^ida + ieAa - m)] (2) 
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FIG. 1. Schematic illustration of Casimir interaction between 
TIs with opposite topological magnetoelectric polarizability 
Oi — —O2. We assume the thickness of magnetic coating is 
much smaller than the distance between TIs. 



where a = 0,x,?/; 7^ = <j^, 7^ = ivpcr^^ 7^ = —ivpcr^ 
and ^/S = ?/^'''7^. a^^y^^ are the three Pauli matrices of the 
spin, and vf is the Fermi velocity of the surface fermion, 
it takes different values for different materials ' , for ex- 
ample, Vf = 1.3 X 10"^ for BisTeg, vf = 1.7 x 10"^ for 
Bi2Se3, in this paper, we take vf = 1.0 x 10~^ for nu- 
merical calculation. Aa are the first three components 
of the electromagnetic potential, m is the surface band 
gap opened by magnetic coating, m = ±\m\ correspond- 
ing to ^ = ±7r. The generalization to = (2n + l)7r is 
straightforward by introducing multi-fermions on TI sur- 
faces. For analytical derivation, we only consider the case 
= ±7r, the general case will be considered in Sec. IV. 

Formally, one can integrate out the contribu- 
tion of surface fermion to get an effective action 
of electromagnetic field on TI surface, Seff{A) = 
— ilndet [i7^(9a + ^eA^) — m]. Up to one- loop approx- 
imation, the quadratic term of effective action Seff{A) 
can be written as: 



s, 



eff 



duj d k 
2n (27r)2 



A(fe,w)n«''(fe,w)Ab(fe,u;). 



(3) 

For the detailed derivation of polarization operator 
tensor 11"'' at finite temperature, we work in Matsubara 
imaginary time formalism: 



E 



(27r)2 



TV[7°G(fc + p)7''G(p)], (4) 



where j3 = l/ksT is the inverse of temperature and 
is the Boltzmann constant, k = {iun^k)^ p = {iuOm^p)^ 
G{p) = i/{j^ka + m) is the propagator of the surface 
fermion. iujri = 2n7rz//3 and iujrn = (2m + l)7rz//3 are the 
finite temperature frequency of electromagnetic field and 
surface fermion respectively. 



The action of surface fermions are relativistic and 
satisfies Lorentz symmetry (if we set the Fermi velocity 
Vf = 1), a similar action and corresponding polarization 
tensor have been considered in graphene system-^^ and 
3-dimensional quantum electromagnetic dynamics , the 
only difference here is that we have only one specie of 
Dirac fermion here, so that the topological parity odd 
term is preserved, after derivation, we find the polariza- 
tion tensor can be divided into three parts: 

U{k) = 2^iUsi{k) + 2v^F^2Us2{k) - mAs{k), (5) 

where Ilsi{k) and n5'2(^) are parity even and IlAs{k) 
are parity odd, their exact forms are: 



iidjiky 




(6) 



where $i, $2 and (p are three parameters, which 
can be derived from Eq.4 straightforward via Feynman 
parametrization and redefining the integration variable 
I = p + xk: 



*1 = -T^^ 



fc2 /? 



dx 



dH 



J (27r)2 

(m^ + vjr {l^ - X{1 - x)k'^) - LOrn (^m + W„)) 



2^,2 



$2 = 



e^v 



-y 



■E 



dx 



XUJnV^vlP^/^Y 

dPl 2x{l — x)v 







(2^)2 [( 



- xou. 



„)2 + 4Z2 + A]'^ 



d^l 



'^''J (2^)^ [(a 



2vpm 



- XLOn) 2 + V%P + A] 



r,(7) 



where A = m^+x(l-x) (cj^ + vj^^ {k^. ^ky)). One can 
carry out the integration over momentum and summation 
over frequency, and get the form of these parameters with 
only the integration over Feynman parameter x\ 



$1 = lo ["^^ tanh(A+) - log (2 cosh (A+) ) 



(A+^A-,/+^r)} 



^2 = 



a I ^ ^ x(l — x) 
dx- 







[tanh(A+) + (A+ ^ A")] , 



ma 



^'^[tanh(A+) + (A+^A-)], 



(8) 



where a is the fine structure constant and 
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r = 



A ± xiujn 



>)/(2r), 



^ _ 2A =F (1 - 2x)iLu„VA - 2x{l - x)v% [kl + kf' 



atVa 



(9) 



By using the series expansion tanht = 1 

2Ej°li(-lVe"'* and log(2cosht) = t-E,°li 
one can rewrite these three parameters as 



(-l)^ ^-2t 



^1 



Jo a/A 4/^2 Jo 



Re/^ — — ) cos(2j7rnx) + Im/^ sin(2j7rnx) 



$2 = 



dx 



x{l — x) 



1 + 2 ^(-l)^e-^'^/^ cos(2jn7rx) 



am 



[ dx—= 

Jo Va 



1 + 2 ^(-l^e-^'^/^ cos(2jn7rx) 



(10) 



I 

where Re/+ and Im/+ take the real and imaginary ihustration in Fig.l), 0(t) is the Heaviside unit step 
part of /+. It is easy to check that in the low temperature function, z = 1, 2 and Zi = 0, ^2 = a. is the value of (j) 
limit (T 0), these expressions coincide with the zero for different surfaces, without lose of generality, we have 
temperature results. By using Eq.3, we get the effective assumed the absolute value of surface band gap equal 



Lagrangian of surface electromagnetic field: 



Stt 



to each other on the two surfaces so that different signs 
of surface band gap corresponding to different signs of 
the topological term aOE • B/iAi^Y in the topological 
quantum field description of TIs. We also note that the 
effect of finite temperature has been implicitly included 
in parameters <l>2 and (j). 



where Fah = daA^ — d^Aa is the electromagnetic field 
tensor. 

Combine Eq.l and Eq.ll, we obtain an action of the 
whole system: 



S ■ 



J dtd^x |-^(^' - B^)0{z)0{a - z) 




- S{z,) 



(12) 



where a is the distance between TIs, and we have 
omitted the thickness of magnetic coating(see schematic 



III. MODIFIED MAXWELL EQUATIONS AND 
CASIMIR INTERACTION 



The Euler-Lagrange equations of the action 12 give the 
Maxwell equations of electromagnetic field with surface 
corrections: 



-V D 



47r 

47r V dt 

5{z - Zi) 



-5{z - Zi) 
-{V X H) 

47r 



X H 

27r dt 27r 



V • B = 0, 
dB 

^ + (vx£;) = o, 



(13) 
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where D = sE and H = B / ji are electric displace- 
ment field and magnetizing field; Ei = CijEjii^j = x^y). 
From these modified Maxwell equations, we get the dis- 
continuous boundary conditions: 

D, {zf) - D, {z-) = -^iB, + 2^1 V • E, 

{zt) - {z-) = ^^E, - 2 (^i^t^^ + ^2VFdxB,) , 

Hy (Z+) - Hy (Z-) = ^,Ey + 2 ( ^ 1 " ^S^^^^.) ^ 

where zf" means ± 0. The other three components 
Ex, Ey and Bz are continuous on the interface. From 
the discontinuous boundary conditions we find that a 
TE mode injection will induce both TE and TM mode 
reflection/refraction. The electromagnetic waves with in- 
jection TE mode in the vacuum can be written as: 

-E' — (1 H~ '^ee)^{ kyCx H~ kxG.y^ -\- Tern{ kzk k G^;)? 
B = {-kzk + k'^Bz) + Teeikzk + k'^Bz) 

y^x H~ kxBy^, (-^^) 

where Tee and rem are reflection coefficients of TE and 
TM mode respectively, the refracted light with refraction 
coefficients tpp and tpm in the TI take the forms: 



E = teeOj{-kyBx + kxBy) + ctem{Pzk - k2Bz), 

B = tee{-Pzk + k'^Bz) + ^^Uj{-kyBx + kxBy), (15) 

where c is the relative velocity of light in TI bulk, k = 
kx^x ~l~ kyBy, k = kj, -\- ky, and Pz is the z component of 
the wave vector in the TI. From the boundary conditions 
we deduced the following equations on the jth boundary: 



/rrrr^em ^ '^z' em ~ 



p,t^£ + ^hirg^ - 1) = -Ai(/.ifc.r(£ 



(16) 



For the TM mode injection, one can write similar equa- 
tions with reflection coefficients rme, ^mm and refraction 
coefficients tme^ trnm- Their solutions are given by (we 
need only the exact form of reflection coefficients): 
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where the superscript j means the jth interface, and 



D={l + s + cP')+(e^ + ^)+ 2$i (fc, + p.) 

\ Pz f^z / 

+2f— + ^+2$iV$iw2^$24fc')- (17) 

\Pz I^Z J 

We note that we have already translated the expres- 
sions into Matsubara imaginary time formalism and we 
assume the influence from permeability can be omit- 
te d, 11 = 1. In imagi nary frequ ency formalism, kz = 
yjuo'^^k'^ and Pz = \J + • For practical calcula- 
tion, we need a form of frequency- dependent dielectric 
permittivity, which can be modeled by'^^''^^ 



K 

J=l 



9J 



(18) 



with K oscillators and for each oscillator, the oscilla- 
tion strength is gj and oscillation frequency is cjj, 7j is 
the corresponding damping parameter. We consider only 
one oscillator and omit the contribution from damping 
parameter here, the generalization to multi-oscillator and 
non-zero damping parameters is straightforward. Then 
the Casimir energy density at finite temperature can be 
deduced from Lifshtz formula: 



E^ 
A 



(19) 

where the prime in the summation means for the n = 
term there contains a prefactor ^, and VS^''^^ are Fresnel 
coefficient matrices on the surfaces, which take the forms: 



l3) 



I em 



I me I mm 



(20) 



IV. RESULTS AND DISCUSSIONS 

It is hard to obtain the full analytical expressions of ^i, 
<I>2 and 0, and a general analyse of the Casimir energy 
for finite surface band gap at finite temperature seems 
to be very difficult. Contrast to the usual calculation of 
Casimir interaction at finite temperature, here the finite 
temperature correction can be divided into two part, the 
one part is the difference between integration and dis- 
crete summation, the other part is from the finite tem- 
perature correction of ^i, <I>2 and (/), and the widely used 
Abel-Plana formula*^^'^^^ does not work here because the 
integration kernel do have singularities on the right-half 
complex plane which have been implicit contained in the 
integral form of <I>i, <I>2 and (j). 
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FIG. 2. (a) Casimir Energy as a function of dimension- 
less distance between TIs at different temperatures, here 
the dimensionless oscillation strength has been chosen to be 
9 — \f9J I^J — 0.45, surface band gap |m| = 500c(;j, and 
T' — T/Tq^ where the definition of To has been given in 
the context, (b) Casimir Energy as a function of dimen- 
sionless distance between TIs for different oscillation strength 
9' — \f9J I^J^ here temperature T = 2To and surface band 
gap \m\— SOOcjj. 



FIG. 3. Equilibrium distance of Casimir force as a function of 
temperature for (a) different surface band gap and (b) differ- 
ent topological magnetoelectric polarizabilities. The dimen- 
sionless oscillation strength has been chosen to be s/gl jooj — 
0.45, definition of dimensionless surface band gap and temper- 
ature are given by = |m|/c<;j and To = ujj /2n respectively. 
Insert shows the detailed difference between m' = 500 and 
infinity surface band gap. 



Here, we are only concerned with the critical surface 
band gap for repulsive Casimir interaction. As shown 
in Ref.^^, the critical surface band gap is much greater 
than room temperature, rric ^ 300i^, so low temper- 
ature expansion is a good approximation. In practical 
calculation, we sum over the first several terms of Eq.l9 
and use the integration over the rest regime to approxi- 
mate the summation with corrections evaluated by Euler- 
Maclaurin formula . 

Before the detailed discussion of results obtained, we 
make a note on units chosen in this paper, we have set 
the Plank constant and velocity of light in vacuum to 
1, and we choose oscillation frequency ujj as the unit of 
energy. The unit of Casimir energy and temperature are 
Eq = cJj/Stt^ and Tq = ujj/27r respectively. 

First, we obtain Casimir energy as a function of dis- 
tance between TIs for different temperature, in Fig. 2 (a). 



and different oscillation strength, in Fig. 2(b). Fig. 2 (a) 
is one of the major results, which shows that increas- 
ing temperature will depress the repulsive peak and re- 
duce the distance between local maximum and minimum 
points of Casimir energy, and at the critical temperature 
Tc, they equal to each other and the repulsive Casimir in- 
teraction vanishes, for given parameters y^/cjj = 0.45 
and \m\ = 500a; j, the critical temperature Tc ~ 2.65To, 
as shown by the red circle dotted line. Such a result is 
well understood because in the high temperature limit, 
Casimir interaction will tend to the classical limit and 
the majority contribution is the zero-frequency term and 
the quantum fluctuation from surface Dirac fermions is 
suppressed. Repulsive Casimir interaction is also quanti- 
tatively influenced by oscillation strength of electromag- 
netic wave in TIs bulk, small oscillation strength will 
decrease the attractive Casimir force from TI bulk and 
profitable for repulsive peak, as shown in Fig. 2(b). 
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FIG. 4. Boundary of repulsive and attractive Casimir inter- 
action as a function of dimensionless distance d = aujj and 
product \m\a for (a) different temperature T /Tq^ (b) different 
oscillation strength g' = ^/gjooj and (c) different topological 
magnetoelectric polarizabilities Q. 

Second we give the local maximum and minimum 
points of Casimir energy(they are both equilibrium dis- 
tances of Casimir force) as a function of temperature for 
different surface band gap, in Fig. 3(a), which shows that 
larger surface band gap will make the critical tempera- 
ture higher and repulsive distance larger. However, such 
a exertion seems to be difficult to achieve and produce 
little effect compared with increasing topological mag- 



netoelectric polarizability. We also give the local maxi- 
mum and minimum points of Casimir energy for topolog- 
ical magnetoelectric polarizability 6> = Stt by introducing 
multi-fermions on TI surfaces, in Fig. 3(b), which shows 
that large topological magnetoelectric polarizability will 
remarkably increase the scope of repulsive Casimir force. 

Then, as a competition, we also give the equilibrium 
distance as a function of temperature for infinite surface 
band gap limit |m| oo, as in Fig. 3(b), which shows 
that, at low temperature, the larger equilibrium distance 
of Casimir interaction for a large but finite surface band 
gap |m| = 500cjj is very close to the equilibrium dis- 
tance of Casimir interaction for infinite surface band gap, 
however, at the critical temperature, Tq ~ 23.5To, the 
larger equilibrium distance and smaller equilibrium dis- 
tance equal, so Casimir force will always attractive when 
T>Tc. 

Finally, similar to the zero temperature case, we also 
give the boundary of attractive Casimir interaction 
and repulsive Casimir interaction as a function of 
dimensionless distance d = acjj and product \m\a for 
different temperature, different oscillation strength and 
different topological magnetoelectric polarizabilities, in 
Fig. 4, and find that the critical product mc a ~ 1/2 is 
still valid, which show that in the short distance limit 
Casimir interaction is dominated by surface fermions, 
where the topological response of surface fermions gives 
a repulsive Casimir interaction and the electromagnetic 
dynamical response of surface fermions give an attractive 
Casimir interaction, these contributions have the same 
magnitude. 



V. CONCLUSIONS 

In this paper, we calculate the Casimir interaction be- 
tween TIs with opposite topological magnetoelectric po- 
larizability and finite surface band gap at finite temper- 
ature, and find that, finite temperature will quantita- 
tively affect Casimir interaction, if Casimir interaction is 
repulsive for proper distance at zero temperature, rising 
temperature will depress the repulsive peak and at a crit- 
ical temperature, the Casimir interaction will be attrac- 
tive for any distance between TIs. We also find that the 
estimation relationship rricCi ~ 1/2 for critical repulsive 
Casimir interaction is valid for different temperature, dif- 
ferent oscillation strength and different topological mag- 
netoelectric polarizabilities, which is useful for practical 
research of repulsive Casimir interaction between TIs. 
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